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Abstract: We study the nuclear symmetry energy of dense matter using holographic 
QCD. To this end, we consider two flavor branes with equal quark masses in a D4/D6/D6 
model. We find that at all densities the symmetry energy monotonically increases. At 



small densities, it exhibits a power law behavior with the density, E, 
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P 
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1. Introduction 

Nuclear symmetry energy is one of key words in nuclear physics as well as in astrophysics. 
Its density dependence is a core quantity of asymmetric nuclear matter which has important 
effects on heavy nuclei and is essential to understand neutron star properties. Although 
much efforts have been given, it is still very poorly understood especially in the supra- 
saturation density regime, see ||, ^ ^ ||, ||, f?J for a review and for a recent discussion. 

From experimental side, the available data do not constrain much the value of the 
symmetry energy at supra-saturation densities. Recently, using the FOPI data on 7r~ /tt + 
ratio in central heavy ion collisions, Xiao et al. j|] obtained a circumstantial evidence for 
a soft nuclear symmetry energy at p > 2p$, where the nuclear symmetry energy increases 
with the density up to the saturation density po and then starts to decrease afterwards. 
Theoretically, almost all possible tools were employed to study the density dependence of 
the symmetry energy. While they showed similar behaviors up to the nuclear saturation 
density, at supra-saturation densities, all possible results one can imagine were predicted 
and no consensus could be reached: some showed stiff dependence (increasing monoton- 
ically with density), while others showed soft one, see Fig. [l| for a typical example. See 
also H for a review. Given this situation, it would be very interesting if we can examine 
the behavior of the nuclear symmetry energy at high densities with a reliable calculational 
tool. 



The gauge/gravity duality [1C, [ll|, |12[ provides a new tool to study strongly interacting 
dense matter, and a few models for QCD [|l4|, based on the duality were constructed. 
Although the true holographic dual of QCD is yet to be constructed, it is worthwhile to 
find out what the new tool says about QCD using available models mimicking the dual 
of QCD. A way to treat the dense matter in confined phase was suggested in and 



a model for transition from nuclear matter to strange matter was proposed in 17]. The 
purpose of this paper is to calculate the symmetry energy of nuclear matter in this model. 
We will find that the symmetry energy is increasing with the total charge Q, showing that 
the symmetry energy of our system has a stiff dependence on the density. Also, we will 
explicitly calculate the density dependence of the symmetry energy at low density to show 

£ S ym ~ P 1/2 - 
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Figure 1: (Color online) Example of density dependence of the nuclear symmetry energy, taken 
from M. Depending on the value of the parameter x, various high density behaviors are possible. 



2. Nuclear symmetry energy in D4/D6/D6 model 

The nuclear symmetry energy is defined as the energy per nucleon required to change isospin 
symmetric nuclear matter to pure neutron matter. In the Bethe-Weizsacker mass formula 
for the nuclear binding energy, it represents the amount of binding energy that a nucleus 
has to lose when the numbers of protons and neutrons are not equal. The semi-empirical 
mass formula based on the liquid drop model has the form: 

E B = a v A - a a {N - Zf/A - a c Z 2 /A 1/3 

-a s A 2 / 3 ±a 5 /AV 4 . (2.1) 

Here Z (N) is the number of protons (neutrons) in a nucleus. The first term is called the 
volume energy since the volume of a nucleus is proportional to A, where A is the total 
nucleon number. The origin of this volume term is the strong nuclear force. The second 
is known as the asymmetry term, which defines the symmetry energy. If there were no 
Coulomb repulsions between protons, we would expect to have equal number of neutrons 
and protons in nuclei in general. The term with a c accounts for the Coulomb interaction 
of all pairs of protons in the nucleus. The last two terms represent the surface energy and 
pairing effect, respectively. Using data for nucleus binding energies, one can determine a 
set of coefficients in Eq. fl2.1|) . 

Due to the invariance of nuclear forces under neutron-proton interchange, iso-scalar 
quantities in a nuclear system are function of only even powers of the asymmetry factor a 
defined by a = (N — Z)/A. Then we can express the energy density per nucleon E(p,a) 
as 

E(p,a) ~E(p,0) +S 2 (p)a 2 , (2.2) 

where p is the nucleon number density and S 2 (p) = 5^^)5=0 is the symmetry energy. 

Now we study the symmetry energy in the D4/D6/D6 model with baryon vertices 
which consist of compact D4 branes and fundamental strings [171- The gluon dynamics 
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is replaced by the gravity sourced by the N c color D4 branes, and two probe D6 branes 
are used to describe the up and down quarks. The bare quark masses are the distances 
between the N c color D4 and two D6's in the absence of the string coupling. 
We can write the metric of the confining D4 background as 

ds 2 = (U/R) 3/2 (-dt 2 + dx 2 + f(U)dx 2 4 ) + (R/U) 3/2 {U/£) 2 (d£ 2 + fdnf) , 

where f{U) = 1 - (U KK /Uf and {U/U K k? /2 = (£ 3/2 + T 3/2 )/2 = £ 3/ V/2. 

We wrap the compact D4 brane on S 4 which is transverse to the color D4 brane. Due 
to the Chern-Simons interaction with RR-field, U(l) gauge field is induced on the D4 brane 
world volume. The source of gauge field is interpreted as the end point of fundamental 
strings. Substituting the equation of motion for gauge field to the Dirac-Born-Infeld action 
of D4 brane with the Chern-Simons, we get the Hamiltonian for the compact D4 brane as 

H m = nj MyJjV\? + ? 2 )yJ ' D{0? + sin 6 d , (2.3) 

where T4 = ^j3^4^3g^ 1 R 3 UKK, D(8) = —2 + 3cos# — cos 3 6 and the prime denotes the 
derivative with respect to 9. We assume that the radial coordinate £ depends only on the 
polar angle 6 of S 4 . 

The fundamental strings out of the compact D4 branes are attached to two D6 branes, 
and they provide the source of U(l) gauge field on the D6 brane. By taking the Legendre 
transformation for the gauge field, we obtain the Hamiltonian which controls the brane 
configuration with fixed charge. 

U m =r 6 J dp^/TT¥^f {q 2 + p 4 4 /3 ) , (2.4) 



where rg = \j-L&V^29^ 1 Uj <:K . Q is dimensionless and related to the number of fundamental 
strings Q by Q = 2 ^% K ^' T6 • Baryons are represented by compact D4 branes, and each of 
them has _/V c fundamental strings attached. Such configuration of compact D brane plus 



fundamental strings are called baryon vertex [13|. The other ends of fundamental strings 
are attached to D6 branes. Therefore, D6 branes are pulled down and compact D4 brane 
is pulled up. As discussed in [16], the length of the fundamental strings becomes zero since 



the tension of the fundamental strings is always larger than that of D-branes. Finally, 
the position of the cusp of D6 branes should be joined to that of the compact D4 brane. 
We consider Q\ fundamental strings attached to one of the D6 branes and Q2 strings to 
another D6 brane. The final configuration is drawn in Fig. (|2|). We denote the slope at the 
cusp of each brane as and yc . The force at the cusp of D6 branes can be calculated 
to give 

dH(Q 2 ) m 



F = 9H(Qi) D6 



dU r . 



dU c 

a 



F [ nl{Qi) + F {2 l{Q 2 ). (2.5) 
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Figure 2: (Color online) Embedding of D-branes with a ^ 0.5. The asymptotic heights of two 
branes are the same (mi = mi = 0.1). The two branes meet at infinity since m\ = m,i. Red curve 
denotes to the position of Ukk ■ 



To make the system stable, following force balancing condition should be satisfied; 

§-F I h = F^(Q 1 ) + F^(Q 2 ), (2.6) 

where Q\ = (1 — a)Q and Q2 = aQ with < a < 1, and Fd& is the force at the cusp due 
to the compact D4 brane. Note that a = (1 — a)/2. 

To find the ground state of our system, we need to consider the energy minimization 
together with the force balancing condition. The total energy of our system is 

E tot = J-^ D4 + H m (Qi) + U D6 (Q 2 ) . (2.7) 

With given Q and a we can calculate the configuration satisfying the force balancing 
condition. The energy density per nucleon E(p, a) is given by Etot in Eq. (2/7) divided by 
the total baryon number Q/N c . 

The symmetry energy is the coefficient of leading term of a and is a function of Q. 
Since Q is proportional to the density of the quark or baryon, this way we can calculate the 
the symmetry energy as function of density. If we further impose minimization of the total 
energy, we can determine the value of a as a function of Q. For m^/mi 7^ 1, there exists 
a transition from a matter with o = 0toa/0ata finite value of Q. This is identified 
as a transition from nuclear to strange matter. For very large Q, a saturates to 0.5, as 
expected. If we take mi = rri2 and if we do not consider isospin violating interactions or 
electromagnetic interactions, then the ground state of the matter would be always with 
a = 1/2. 

The explicit form of symmetry energy per nucleon can be written as 



N B J (Q2 + 4w 8/ 3 p 4)3/2 ' 

where y is the embedding solution of D6 brane with a = 1/2. Notice that Nb = Q/N c 
and so the symmetry energy ( |2.§| ) contains N c factor. We need to factor this N c out for 
the reason we discuss in Section ||. Our results are given in Fig. ||| 
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Figure 3: (Color online) Solid line is the calculated symmetry energy as a function of the density. 
For illustration purpose, we take A = 6 and Mkk = 1-04 GeV. The dotted line is for S 2 ~ p 1 / 2 . 



Note that so far we use p for both the coordinate and the density. Hereafter p is only 
for the density. Although our main point here is the calculational scheme of the symmetry 
energy rather than the exact numerical numbers, it is better to see how the numbers fit to 
the reality. To fix the energy scale, we rewrite the energy density per nucleon as 

where E = E to t/Te. One may determine the values of A and Mkk by using the quark 
mass and meson mass as inputs into two basic relations: m q = XMkkVooI / (2 2//3 9vt) and 
AfJ = 0A6M kkVoo with our coordinate choice. From the the non-anomalous rj' mass, 
~ 390 MeV pofl , and the quark mass, M q ~ 41 MeV, 1 we can determine the parameters 
of the D4/D6 model: M KK = 1.04 GeV (2§ and A = 6. In this case, Q ~ 1.2 corresponds 
to the normal nuclear matter density pq. 

Now, we lay emphasis on two aspects of our results, which are rather insensitive to the 
choice of A and Mkk- One is the stiffness of the symmetry energy S2 in supra-saturation 
density regime, and another is its low density power law behavior S2 ~ p 1 ^ 2 . 

The power law behavior of S2 m low density can be understood by calculating ana- 
lytically in a special limit, m q -4 00 and p — > 0. In this case, the solution of D6 brane 
embedding becomes trivial, y = and we can integrate ( |2.8|) analytically to have 



The current experimental result of the symmetry energy can be summarized by a fitting 
formula 

s 2 ( P ) = c(p/p y (2.11) 



1 As well-known, the quark mass in D4/D6 model could be different from that in QCD by a constant 
factor. To obtain the constant we need to compare the scalar two point function obtained in D4/D6 model 
with that in the operator product expansion of QCD, for example see Jjlj ]. 
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with c ~ 31.6 MeV and 7 = 0.5 — 0.7 in the low density regime, 0.3po < p < pcb see || 
for example. With our choice of A, Mkk, we obtain 7 ~ 0.5 and c ~ 27.7 MeV from Eqs. 
J2.8|) and ( |2.9| ) . Notice that the value of 7 in our results is rather insensitive to the value 
of A and Mkk, while the value of c depends on them. 

One can understand the stiffness based on the property of the branes: suppose two D6 
branes meet with compact D4 at a point with different polar angles 6\ and #2 as drawn in 
Fig. ||. Here Q\ and Q2 with Q\ + Q2 = Q fixed are the number of the strings attached to 
each D6 brane. The force balance condition is nothing but the minimization of the action 
with respect to the position of the contact point. To balance the pulling force of the compact 
D4, the position of the contact point is located when each D6 brane balances roughly half 
of the pulling force, which is a statement supported from the numerical analysis. Since the 
upward force is proportional to the product of effective tension, ~ tqQ, times the projection 
factor to the vertical direction, cos 9, we have Qicos#i ~ Q2COS#2- If #1 > #2> Qi > Q2- 
Namely, more strings should be attached to the lower brane in the figure 2. Therefore, the 
asymmetry in the number of the attached string is due to the angle difference. Notice that 
each end point of the string provides the 'electric' flux which contributes to the energy 
of the brane. Since the flux of charge Qi is confined in each brane, the total energy is 
f(Qi) + f{Q2)-, where f(Q) = T-Ld&{Q)- Since / is a monotonically increasing function, the 
minimization of the energy with respect to the variation of Q\ requests Qi = Q2 = Q/2. 
As Q increases, the effective brane tension increases and so maintaining the angle difference 
costs more and more energy. Furthermore, since f"(Q) ~ \ j\J~Q is positive, the symmetry 
energy, Q 2 f' '(Q /2) /Nb, is a increasing function. 

The Coulomb repulsion discussed above is, of course, not the electromagnetic one. The 
local U(l) charge is holographic dual to the the global baryon number of the boundary 
theory. However, the repulsion in the dual bulk theory means repulsion in 4 dimension as 
well. From the boundary theory point of view, such repulsion is simply due to the presence 
of the baryon charge. So the origin of the repulsive nature is mysterious from the boundary 
point of view. To understand this, we notice two facts: The first one is that the charge 
carriers are fermions since the charge is introduced by the D4/D6 fundamental string end 
points, not from the bulk R-charge of type IIA gravity. The second is that in the boundary 
theory the tendency oi N = Z by the Pauli principle, while in the holographic dual bulk 
theory it is the Coulomb interaction that requests N = Z. Since two origins should be the 
same, we may suggest the Coulomb repulsion as the holographic Pauli principle. See [18] 
and [ fil| for similar observations in a different context. 

Finally, we study the effect of small isospin violation by considering mi 7^ 7712- We 
find in this case that the symmetry energy is almost the same with the case with isospin 
invariance for the mass ratios of order one. 



3. Summary and Discussion 

In summary, we calculated the symmetry energy of dense matter in the D4/D6/D6 model. 
To obtain the symmetry energy in nuclear matter with charge symmetry, we considered 
the case with mi = rri2 and found that the symmetry energy is increasing with the total 
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charge Q, showing a stiff nuclear symmetry energy. It is universal in the sense that the 
result is independent of the value of A and Mkk- We also studied the low density behavior 
with power 7 to be ~ 1/2, which is again independent of the value of A and Mkk and is 
close to the value suggested by experiments, 7 = 0.5 — 0.7. 

One subtle point we mentioned in the main text was about the factor N c . The reason 
we divided this factor out is as follows. In our model, the same flavor quarks form a nucleon; 
for instance, proton in our model consists of N c mi quarks and neutron has N c rri2 quarks 
in it. Hence, the total number difference of quarks is N c times the number difference of 
neutrons and protons, resulting in the overall iV c factor in the symmetry energy. However, 
in reality, where N c = 3, proton consists of two up quarks and one down quark, and neutron 
contains one up quark and two down quarks. The total number difference of quarks is the 
same as the number difference of neutrons and protons. Therefore, in order to compare 
our result with the realistic case, we have to divide the symmetry energy (|2.8|) by N c , 

We also studied the effect of isospin violation by considering mi ^ 771,2. However, the 
symmetry energy in this case turned out to be almost the same with the case with isospin 
invariance for the mass ratios of order one. 

Now, we list some generic cautionary remarks. In holographic approaches, it is not 
clear how to encode attractive scalar contributions that are essential to describe nuclear 
matter. In addition, for extreme large density limit, the back-reaction effect from the dense 
matter is not negligible, and so it should modify our result. 

Finally, we comment on a future investigation. In conventional approaches, for instance 
sec ||, ^5| , at very low densities p <^ po, the dominant contribution to the symmetry energy 
is coming from the kinetic energy which encodes the Pauli principle. This is because the 
kinetic contribution to the symmetry energy is ~ p 2 ^ 3 , while the one from interactions 
starts from ~ p 1 due to the linear density approximation which works well at very low 
density. The origin of the factor 7 = 2/3 is the dispersion relation E ~ p 2 together with 
the sharp Fermi surface. In our case, the fact 7 = 1/2 suggests that either the dispersion 
relation is anomalous like E ~ p 3 / 2 or Fermi surface is fuzzy |24j due to the strong 
interaction. This poses an interesting future study. 
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